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An electrical method is used to study the early stages of coalescence of two low-viscosity drops. 
A drop of aqueous NaCl solution is suspended in air above a second drop of the same solution 
which is grown until the drops touch. At that point a rapidly widening bridge forms between 
them. By measuring the resistance and capacitance of the system during this coalescence event, 
one can obtain information about the time dependence of the characteristic bridge radius and its 
characteristic height. At early times, a new asymptotic regime is observed that is inconsistent with 
previous theoretical predictions. The measurements at several drop radii and approach velocities 
are consistent with a model in which the two liquids coalesce with a slightly deformed interface. 

PACS numbers: 47.55.D-, 68.03.-g, 47.55.df, 47.55. N- 



I. INTRODUCTION 

All around us we see fluid drops joining together: rain- 
drops splash into a pond and fuse with it; individual 
drops falling from a faucet merge together to fill a glass 
of water. It is easy to forget the wonder of such a ubiq- 
uitous phenomenon. There is a change in topology when 
two fluids coalesce. As soon as they come into contact, a 
fluid bridge is formed between the two masses. The ini- 
tial radius of the bridge is much smaller than the macro- 
scopic dimensions of the flow. Interfacial tension then 
widens it until the two drops merge into a single en- 
tity. A video sequence of this process is shown in Fig. 
1. Another common example of a topological transition 
is the inverse process to coalescence, that is, drop break 
up[TJ |5J [21 miS]. There, a single mass of fluid separates 
into two segments joined by a thin neck. In that case 
the topological transformation proceeds as a physical di- 
mension, the neck radius, approaches zero causing the 
dynamics to approach a singularity. When two drops 
coalesce, we expect similar singular behavior. 

Such fluid transitions have often been compared to crit- 
ical thermodynamic phase transitions, as this separation 
of length scales often leads to universal behavior [6l III [8]. 
Although it is an appealing and useful framework, it was 
recently discovered that not all fluid-breakup singulari- 
ties obey universal dynamics jH [10]. In light of this, it 
is imperative to consider other familiar fluid transitions, 
such as drop coalescence, to see if they, too, behave in un- 
expected ways. Moreover, drop coalescence is of practical 
as well as purely scientiflc importance. Viscous sintering, 
emulsion stability and mixing in microfluidics often need 
to be controlled in industrial processes. In this paper, I 
employ an electrical method to explore the drop coales- 
cence transition at low viscosities at times three orders 
of magnitude earlier than previous optical experiments. 

Coalescence processes occur in both the viscous regime, 
where the primary force opposing the widening of the 
bridge between the drops is due to viscous dissipation, 
and the inviscid regime, where the widening of the bridge 
is opposed primarily by inertial forces. In the inviscid 
regime studied here, the radius of the fluid bridge be- 




FIG. 1: Two drops of aqueous NaCl solution at saturation of 
radius A = 1 mm are coalescing. The frames are 69 /xs apart. 
The white spot in the bridge center is an optical artifact due 
to the drop lensing the light source located behind it. 



tween the two coalescing drops, r, is much greater than 
the viscous length scale of the system, — /i^/p7, where 
fi is the dynamic viscosity of the fluid, p is its density, 
and 7 is the surface tension. For water coalescing in 
air, « 14 nm, and viscous effects can be neglected for 
much of the coalescence. This regime has been studied 
less than its high- viscosity counterpart [11] [T5] because 
the very rapid initial motion of the low viscosity fluid 
is difficult to resolve in experiments and computations. 
However, theory has predictions. 

A straightforward scaling argument |llj can be used 
to describe coalescence in the inviscid regime. To ini- 
tiate coalescence, the drops must be brought very close 
together. Soon after the bridge is formed, the gap width 
between the two drops, d, will satisfy d <ti r. In this case, 
a balance between surface tension and inertia leads to 



pa dt 



(1) 



It is assumed that if the two droplets are brought to- 
gether sufficiently slowly, then they will maintain a spher- 
ical shape. For hemispherical drops, d — r'^/A, where A 
is the drop radius, as shown in Fig. 2. The resulting dif- 
ferential equation can be solved, where to is the instant 
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at which coalescence occurs, and c is a proportionality 
constant of order unity: 



(2) 



(a) 



This scaling law has been supported by simulations 
studying the coalescence of low-viscosity fluid drops in 
the absence of an outer fluid[T3J[T3]. However, due to the 
speed and the geometry of the transition, experimental 
studies have been unable to confirm the applicability of 
this scaling law to times {t — Iq) = t < \Q fis. 

Previous experiments have observed r oc (t — to)^^^ for 
(t — ^o) > 10 MS, using high-speed imaging at rates up 
to 10® frames per second [121 US] [T7] , as well as ultrafast 
x-ray phase contrast imaging[18' . However, one aspect of 
the data suggests that the dynamics may not behave as 
we would expect. In one experiment, a small DC voltage 
was placed across the drops, and it was found that the 
initiation of electrical contact occurred 20 — 50 fis before 
coalescence could be observed visually [15]. An electri- 
cal method introduced by Case and Nagel [TH], expands 
the measurement range down to r ~ 10 ns, a region cur- 
rently inaccessible to imaging experiments. Results from 
this data closer to the instant of coalescence indicated a 
new asymptotic regime not predicted by the scaling ar- 
gument given above. Moreover it suggested a solution to 
the discrepancy found between the instant of electrical 
contact and the apparent initiation of coalescence. 

Electrical methods for studying drop coalescence have 
been used in other experiments. The electrical method 
used in Case and Nagel's experiment is similar to one 
previously developed by Burton et al[5U], which used a 
small DC voltage to measure the resistance of a mercury 
droplet during break up. Case and Nagel extended Bur- 
ton et al's technique by using an AC voltage to measure 
separately both the time-dependent resistance and the 
capacitance of two coalescing drops of aqueous NaCl so- 
lution. This enabled them to infer the geometry of the 
coalescing region as early as 10 ns after the instant of co- 
alescence. They found a new asymptotic regime at early 
times that is not consistent with the predictions of the 
simple scaling argument outlined above. This behavior 
occurs for r < 10 /xs so that it is entirely in the region 
that cannot be studied by direct imaging. In addition, an 
AC electrical method was used by Lukyanets and Kaveh- 
pour to study the rest time of coalescing drops [21J . These 
results (at voltage magnitudes three orders of magnitude 
larger than the largest used here) suggest that deforma- 
tions resulting from high electric fields in the gap between 
the two drops may introduce errors in measurements. In 
the experiments described here, varying the voltage and 
frequency by several orders of magnitude does not sig- 
nificantly affect the results. This is discussed further in 
appendix A. 

In this paper, I expand upon these measurements and 
provide a more detailed experimental description. I vary 
experimental parameters such as the drop diameter and 
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FIG. 2: (a) Experimental Setup. Two acrylic tubes of length 
5 cm long and inner diameter 0.95 cm are secured in line with 
one another. Changeable nozzles of radius A are attached to 
each tube. The nozzle tips are separated by 2A. A drop of 
aqueous NaCl solution is formed on the upper nozzle using a 
microliter syringe, and the lower drop is then slowly grown 
until the two drops coalesce, using a variable speed syringe 
pump (Kazel R99-FM) with syringe sizes varying from 50 
to 20 mL, and injection speeds varying from 0.21 to to 70.0 
ml/hour, (b) Coalescence of two drops. Two drops of radius 
A meet at a single point. A bridge of radius r and height 
d forms and expands due to the interfacial tension 7. For 
hemispherical drops, d ~ /A. 



approach velocity in order to explore further the surpris- 
ing behavior seen in the initial experiments. My experi- 
ments support the hypothesized new asymptotic regime 
at the earliest times measured. 



II. EXPERIMENTAL DESCRIPTION 

A. Impedance Measurement 

In this experiment, two drops of aqueous sodium chlo- 
ride solution at saturation coalesce in air at room tem- 
perature. At saturation, or 26 % NaCl by mass, the fiuid 
parameters of salt water are: fiuid density p = 1.1972 
g/cm, kinematic viscosity v = 1.662 cSt, surface ten- 
sion 7 = 82.55 dyne/cm, and conductivity cr = 0.225 (fi- 
cm)-i[23]. 

As shown in Fig. 2, two acrylic tubes filled with salt 
water were aligned vertically. Tefion nozzles of radius 
A were attached to each tube, facing each other. Gold 
electrodes were immersed in the salt water at the end of 
each tube opposite the nozzle. A known quantity of fiuid 
was injected into each nozzle, forming two approximately 
hemispherical drops separated by a small distance. An 
AC voltage of frequency / and magnitude \ V\ was applied 
across the electrodes, and the lower drop was then slowly 
grown at a fixed rate until the two drops coalesced. The 
complex impedance of the experimental cell, Zceii, was 
measured as a function of time during the coalescence. 
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FIG. 3: Measurement circuit. Gold electrodes 0.5 mm in di- 
ameter and 1 cm long are secured in the measurement cell 
and connected to the circuit shown. The AC source is an HP 
3325 A function generator (Hewlett-Packard). The upper left 
branch consists of known circuit elements {Rt and Ct), while 
each lower branch is the input impedance of the oscilloscope 
(i?o) in parallel with the capacitance of the coaxial cables 
connecting the circuit to the oscilloscope (Co). In place of 
a traditional oscilloscope we use an NI PCI-5105 high-speed 
simultaneous sampling digitizer. The impedance of the cell 
can be separated into three contributions, added in series: 
Zeiectrodes, Z fluid, and ZcR- Zeiectrode cau be modeled as a 
frequency-dependent capacitance due to the double layer, in 
parallel with an equivalent resistance due to charge transfer, 
as shown in the dotted circle [22]. The impedance of the co- 
alescing region, Zcr can be modeled as a capacitor, due to 
the large exposed surface of the two drops, in parallel with a 
resistance, as shown in the dashed square. 



To evaluate Zceih we used the Whcatstone bridge ar- 
rangement shown in Fig. 3. A known impedance, Zt, was 
connected in series with our measuring device, a National 
Instruments PCI-5105 simultaneous sampling digitizer. 
The effective input impedance of the PCI-5105, Zq, is its 
stated input impedance, Rq, in parallel with the cable ca- 
pacitance Co of the coaxial cable. The experimental cell 
was also connected in series with the measuring device, 
and in parallel with the combination of Zt and Zq. A 
Hewlett-Packard HP3325A function generator was con- 
nected in series with Zf. 

The voltage Vr was measured between Zt and Zq, and 
Vs was measured between Zceii and Zq. These voltages 
were sampled simultaneously at a maximum rate of 60 
MHz. The sampling rate was 10/, where / is the fre- 
quency of the input sine wave from the function gener- 
ator, except for data taken at / = 10 MHz, which was 
sampled at the maximum rate of 60 MHz. The volt- 
ages were read into Labview (National Instruments) and 
analyzed. The analysis averaged the incoming signals 



over a single period to find the ratio of their amplitudes, 
as a function of time, t. In addition, the anal- 
ysis compared the input signals to a known sine wave 
and found the relative phase shift of each signal versus 
t. These phase shifts were subtracted to find, Ac/), the 
phase shift between Vr and V, as a function of t. 
Using the complex equation: 

\Vs\ \Zt + Z,\ ' 

allowed Z^eii to be calculated as a function of the known 
circuit elements Zt and Zq and the measured values of 
\Vr\l\Vs\ and A0: 

2\Vr\ 

ReiZceii) = -^{{Re{Zo) + ReiZt))cosA<l) 
\Vs\ 

-(/to(Zo) + Im{Zt)) sin A(j)) - Re{Za) 
2\Vr\ 

Im{Z,M) = ^^(((i?e(Zo) + i?e(Zt))sinA0 
\Vs\ 

-f (/m(Zo) + Im{Zt)) cos A(j)) - Im{Zo). 

The circuit was calibrated by replacing Zceii with 
known circuit elements, and the measured values of 
Re(Z) and Im{Z) were shown to be consistent across 
the frequency range with the values of the known circuit 
elements. Eq. 3 assumes that the input impedance for 
the PCI-5105 is identical for both input channels. This is 
not necessarily the case, and the analysis allowed this to 
be varied in order to calibrate the cell. Within the known 
error of the input impedances, however, they were iden- 
tical, and this equation is accurate. 



B. Isolating the impedance of the coalescing region. 

The impedance of the experimental cell, Zceiu has 
three distinct contributions, as shown in Fig. 3. It can be 
shown that for the voltage across the cell \ Vceii \ ^ bQmV , 
the interaction of the electrodes with the solution pro- 
duce an effective contribution to the impedance that 
is equivalent to a resistor in parallel with a frequency- 
dependent capacitor [22 . We represent this contribution 
as Z electrodes- Also, the "Coalescing Region" (defined 
as the region between the tips of the two nozzles) con- 
tributes an impedance Zcr- Finally, the fluid between 
the electrodes and the coalescing region contributes an 
impedance Z fluid- Z electrodes depends on the frequency 
/ while Z fluid is independent of /. As these contributions 
are in series, we can write 

Zccll — Zelectrodes + Zfi^id + ZcR- (4) 

If the two nozzle tips are brought into contact, 
Zciosed = Zelectrodes + Z fluid- A representative measure- 
ment of Zciosed as a function of frequency is shown in 
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FIG. 4: Cell impedance. The real and imaginary parts of 
Zciosed = Zeiectrodes + Z fluid are shown as functions of fre- 
quency. The open symbols show Im{Zcioaed), and the closed 
symbols show Re{Zciosed)- 



Fig. 4 for j4 = 1 mm. This measurement allows Zqr to 
be isolated from the other contributions in the cell: 



ZcR — Zcell — Zciosed- 



(5) 



During coalescence, Zcr can be considered as a resis- 
tor, RcR, representing the resistance of the two drops 
and the bridge between them, in parallel with a capaci- 
tance CcR, representing the capacitance of the conduct- 
ing surfaces. Both the resistance and the capacitance 
change with time, and are related to the geometry of the 
coalescing region. 

During the 10 /is before coalescence occurs, Zcr = 
—i/2'KfCcR represents the capacitance of the drop tips 
as well as the other conducting surfaces of the cell. We 
consider Cj„it, the capacitance of the two drops just be- 
fore coalescence occurs, to be in parallel with the capac- 
itanc;e of the rest of the system, CceM- This approxima- 
tion is supported by electrostatic simulations (see Ap- 
pendix B). Thus, in order to isolate Cj„jt, we calculate 
Cinit = CcR. — Cceii, whcrc Cceii IS the measured capaci- 
tance of the cell when no drops have been formed on the 
nozzle tips. In our experiments, we measure both Rcr 
and CcR as functions of time after the bridge is formed, 
as well as Cmit before the bridge is formed. 



III. RESULTS AND DISCUSSION 

A. Resistance and Capacitance during Coalescence 

The resistance of the coalescing region is shown versus 
r = f — to in Fig. 5(a) for A = I mm. We determine the 
instant of coalescence (to) from the phase shift between 
\Vr\ and \Vs\ to within a period of the oscillation, ±1//. 

We adjust tg within this range such that the earliest data 
taken at a given frequency overlaps the data taken at 




FIG. 5: Measured resistance during droplet coalescence, (a) 
-Rcflversus t = {t — to). The solid line shows Rcr = 
1.2 ■ 10~V-^ + 0.8r"i/^ + 23.3. (b) Rcr - 23.3 Q. versus 



r. The solid line shows Rcr = 1.2 



The dashed line shows Rcr 



-1/2 



In each case, A = 1 



mm, and the drops approach one another at a rate of 0.0004 
A/ms. The data is an average of 24 individual coalescence 

events, six obtained at each of four measurement frequencies. 
The error bars reflect the spread in these measurements as 
well as systematic error due to inaccuracies in the measure- 
ment of Z electrode and due to the choice of to, the instant of 
coalescence. Error bars are shown both for r and for Rcr- 



higher /. This is impossible for the highest frequency 
data. 

The data is well described by the form Rcr = aT~^ ■+■ 
(3T~^f'^ -|- 5. The best fit to the data, shown by the solid 
line, gives a = (1.2 ± 0.3) • 10"^, /3 = 0.8 ± 0.2, and 
8 = 23.3 ± 15. Fig. 5(b) shows the same data as in Fig. 
5(a), plotted as Rcr — 23.30 versus r. The dashed line 
shows a power law r~^/^ while the solid line shows the 
same fit as in Fig. 5(a). 

Cinit versus r is shown in the filled symbols in Fig. 
6, and is constant within error. The solid line shows 
the average value of Ccr = 1-30 ± 0.14 pF. The open 
symbols show Cinit = CcR-Cceii versus r, where Cceii = 
0.89±0.02 pF. The average value Ccr = 0.41 ±0.14 pF is 
shown by the dashed line. All capacitance measurements 
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FIG. 6: Measured capacitance before droplet coalescence. 
The closed symbols show Ccr versus r. The solid line shows 
CcR = 1.3 pF. The open symbols show Cinit = Ccr — CcM 
pF. The dashed line shows Ci„it = 0.41 pF. The data is an 
average of 3 individual coalescence events taken at / = 10 
MHz. The error bars reflect the spread in these measure- 
ments as well as systematic error. 
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FIG. 7: High Speed Imaging, r versus r is shown for a repre- 
sentative coalescence event. The solid line shows r = 3.2t'''^°. 
The frame rate is 144,000 frames per second. 



are obtained at 10 MHz, as at lower frequencies, \Vs\ is 
comparable to electrical noise. 



B. High-speed imaging data 

In addition to our electrical measurements, we verified 
previous measurements of the bridge radius during drop 
coalescence using a high-speed digital camera (Phantom 
v. 7) running at 144,000 frames per second. Images are 

shown in Fig. 1. The resolution used was 26 /im/pixcl, 
and we used simultaneous electrical measurements to de- 
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FIG. 8: Two geometries for coalescence, (a) Two hemispher- 
ical drops of radius A coalesce. We separate the resistance 
RcR into three parts. Rhemi is the resistance of the hemi- 
spherical shapes, which are cut off as they come into contact 
with the bridge. Rtridge gives the resistance of the bridge with 
radius r = rtr- (b) Two drops of radius A coalesce with flat- 
tened tips. The radius of the flattened region is given by rfiat- 
For r < rfiat, d = constant, while for r > r/jat, d = /A. 
The resistance of the flattened hemispheres is given by Rdiat- 



termine to. A sample measurement of r versus r for 
A = \ mm is shown in Fig. 7. The best fit to the 
data for r < 1 ms gives r = (3.2 ± o.5)r" '™±o o2 
T > 100/Lts , which is consistent with previous measure- 
ments. This exponent is consistent with the scaling ar- 
gument assuming d oc summarized by Eq. (2), which 
predicts r = 2.3t^/^. 



C. Comparison to predictions from scaling 
argument 

We predict Rcr by considering the geometry of the 
coalescing region. Rcr can be separated into three pieces 
connected in series: Rcr = R 



Rbridge + Ri 

Rupper is the resistance! of the upper drop, Riower is the 
resistance of the lower drop, and Rbridge is the resistance 
of the bridge between the two drops. 

We assume that both drops arc hemispheres of radius 
A, and so Rupper = Riower = Rhemi- Therefore, we write 
Rcr = '^Rhemi+ Rbridge- Duc to the vertical alignment of 
the nozzles, gravity can distort the drop shapes slightly. 
However, for length scales much smaller than the cap- 



upper 
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illary length Ic = \/"f/pg, the surface tension pressure 
maintaining a spherical shape is stronger than gravity. 
For the aqueous NaCl solution, Ic — 2.7 mm, while our 
largest drops have A = 2 mm. 

We model Rtemi as a hemisphere truncated by a plane 
parallel to the flat surface of the hemisphere. The plane 
intersects the hemisphere with a radius of r^^ , as seen in 
Fig. 8(a). We numerically calculate the resistance of this 
shape using the electrostatics calculation package EStat 
(FieldCo), varying rtr over several orders of magnitude. 
We find that Rhemi — ^/^f'tr'^ ; where a is the conduc- 
tivity of the fluid. 

Rbridge Can be estimated directly. The resistance of 
a roughly cylindrical object scales as the length divided 
by the area, yielding Rbridge ~ djai^r^. For two hemi- 
spherical drops, the gap width d ~ /A, and therefore 
Rbridge ^ (o"""^)"^; & coustaut. Using thcsc results and 
the relation between r and r from the scaling argument 
in Eq. (2), we find 



1 



1 



ottA 



(6) 



This equation is calculated for the salt water and air 
system with A = I mm, Rcr = 0.97r~^/^ -I- 14. For 
T > 10 /iS, this prediction is in qualitative agreement 
with our data, as can be seen in Fig. 5(b). However, for 
T < 10 ^s, our data shows Rcr to be significantly larger 
than predicted. Thus the early-time data suggests a new 
asymptotic regime not included in the scaling prediction 
ofEq. (2). 

We can estimate Cinit in this geometry by modeling 
the system as two conducting hemispheres separated by a 
distance z. When z <C A, the capacitance of this arrange- 
ment of conductors is comparable to that of a sphere of 
radius A/2 suspended with its tip a distance z above an 
infinite conducting plane, which can be solved analyti- 
cally. An approximation for z/A<^ \ [MJ shows 



7reoA[ln( 



2z> 



1. 



(7) 



Due to the logarithmic dependence of Cinn on z, un- 
certainty in the measurement of Cinit = 0.41 ± 0.14 pF 
leads to enormous variation in the calculated value of z: 
from 190 nm for C,„,t = 0.27 pF to 8.05 • 10"^ nm for 
Cmt = 0.55 pF. 



D. An Alternative Interpretation: Flattened drop 

tips 

To explain the discrepancy at small r between the pre- 
dictions of the model and the data, a modification to 
the coalescence geometry was proposed [19]. In deriving 



r^/^, it was assumed that d 



However, if the 



drop tips are slightly flattened, as in Fig. 8(b), a different 
dependence is found for r(T) at early times. 



For hemispherical drops with a flattened tip of radius 
Tfiat, d cc only when r > rfiat- For r < rjut, d is 
constant, and the problem is equivalent to that of a hole 
opening in a thin fllm due to interfacial tension. As long 
as we remain in the inviscid regime, Eq. (1) still applies, 
and solving for d constant, we flnd 



(8) 



In this geometry, Rcr is calculated by replacing Rhemi 
for the undistorted case with Rdist, the resistance of the 
distorted hemisphere in Fig. 8(b) that has a flattened 
tip of radius Vfiat- The flattened hemisphere has a small 
opening in its base of radius r which corresponds to the 
bridge. A numerical solution shows Rdist = l/4r(T. We 
estimate Rbridge = djcn^r^. Combining these contribu- 
tions with the time dependence seen in Eq. (8) yields 



2f7 7 T 



cr7r7 



(9) 



Thus, when r < rfiat, RcR is independent of A. 

The predictions of this model can be compared to the 
data shown in Fig. 5(a). A transition time tt from the 
T--1/2 behavior to the behavior is determined to be 
0.87 fis <tt < 6.5 ^s with the best fit being tt = 2.4 ^s. 
For t < tt, Rcr ~ (1-2 • lO'^ ± 3 • 10-'')r-i. From the 
argument above, this prefactor is {pd/^Y^'^ /2a. Com- 
paring the prediction to the data yields d — 200 ±100 
nm. 

The contribution from the r"^ term is negligible. Rcr 
crosses over from r^^ behavior to behavior at a time 

= (0.077 s/cm-3/2)rf3/2^ When d = 200 ± 100 nm, 

= 7 ± 5 ns. This is beyond the measurement window; 
the experiments would only resolve Rcr. oc d^^"^ /t at the 
earliest times measured. 

In this model, Cinit is dominated by the fiattened re- 
gion. We approximate Cinit as a parallel plate capacitor 
of area i^r'^fiat and separation d. Using Equation (2) and 
the crossover time tt from r^^/^ to r^^ behavior gives 
22 /im < rfiat < 59 /xm with the best fit value (from the 
best fit for tt above) rjiat ~ 36 /xm. This approximation 
then yields a capacitance 0.04 pF < C < 0.97 pF with 
the best fit C = 0.18 pF. Using electrostatic simulations 
in combination with Eq. (7), I find that the hemispher- 
ical region contributes approximately 0.07 pF to the ca- 
pacitance. Assuming these parallel contributions can be 
added, I find that 0.11 pF < Ccr < 1.04 pF, which is 
consistent with the measurement Cinit = 0.41 pF. 



E. Varying the drop radius A 

The model with the flattened tips has two regimes. Eq. 
9 should hold for t <^tc, where tc represents the time at 
which Rcr. crosses over from to behavior. At 
longer times, t ^ tc, RcR^ Eq. 6 should apply. Varying 
the drop radius A should only affect data for f:^ tc- 
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FIG. 9: Varying drop radius. Tiic open symbols show A = 2 
mm. The closed symbols show A — 0.75 mm. (a) Rcr versus 
T. Data shown is the average of 12 individual coalescence 
events for each A. Three coalescence events were measured 
with each of four different /, from / = 10 kHz to 10 MHz. 
The solid line shows Rcr = 1.4 ■ lO^'V^^ + 0.9t-"--'° + 40. 
The dashed line shows Rcr = 1.2 ■ lO^'V"^ + 0.7r"" "" + 10. 
(b) Cinit versus r. Data shown is the average of 3 coalescence 
events taken at / = 10 MHz. The dashed line shows Cinit = 
0.52 pF. The solid line shows dnit = 0.73 pF. 



I have measured Rcr for ^ = 2 mm and A = 0.75 
mm in addition to the ^ = 1 mm measurements already 
shown. For a drop with a flat tip, we expect that for 
t > tt, Rcr ~ + 7 for A = 2 mm, and Rcr w 

1.0r~^/^-|-19 for A = 0.75 mm. At early times, we expect 
no change outside of error from the A = 1 mm data for 
both drops. 

The measured Rcr versus r for A = 2 mm and A = 
0.75 mm is shown in Fig. 9(a). For A = 2 mm (open sym- 
bols), the best fit to the data is the dashed line: Rcr. — 
(1.2±0.3) • lO-V-i + (0.7±0.2)t-°-50 + (10± 10). For A 
= 0.75 mm (closed symbols), the best fit is the solid line: 
Rcr = (1.4±0.3)-10--V-i + (0.9±0.2)t-o-5" + (40±15). 
Our measurements thus give prefactors that are qualita- 



tively consistent with those predicted by the flattened tip 
model. 

Cinit versus r is shown for A = 2 mm and A = 0.75 mm 
in Fig. 9(b). Cceii for A = 2 mm is measured to be 0.92 
pF, and Cceii for A = 0.75 mm is 0.81 pF. The measured 
capacitance for A = 0.75 mm (shown in closed symbols) 
is larger than that seen for A = 1 mm, and the average 
value before coalescence is Cinit — 0.84 = 0.73 ± 0.2 pF 
(shown by the solid line.) For A = 2 mm, we observe 
Cinit ^0.92 = 0.52 ±0.14 pF, which is within error of the 
value observed for A = 1 mm. 

In summary, when A is varied, we observe behavior 
which is consistent with the coalescence of two slightly 
flattened drops. For r < 1 /US, no difference is observed 
in Rcr when A is increased by a factor of 2.7. For r 3> 1 
/US, the observations are consistent with the prediction 
that Rcr should increase for smaller drops. In the ca- 
pacitance measurements, we observe increases over the 
1 mm measurement for both A = 0.75 mm and A = 2 
mm. An increase in deformability could account for the 
increase observed for A = 2 mm. For A = 0.75 mm, 
we are unable to reach the low approach velocities used 
for larger drops, which may be responsible for the large 
capacitance observed. 



F. Varying the Drop Velocity 

The drops are brought together at a non-negligible 
approach velocity v. As v is increased, air effects will 
become more marked, particularly for larger drops, and 
may change the drop shape. We isolate such effects by 
varying the approach velocity. To account for the effects 
of different A, our units of velocity are A/ms. 

We show Rcr versus r in Fig. 10(a), where the ap- 
proach velocity is varied by a factor of 5 for A = I mm. 
The closed symbols represent v = 0.0004 A/ms, while 
the open symbols represent v = 0.002 A/ms. Varying 
the velocity by this amount does not appreciably change 
the average data. 

Fig. 10(b) shows Rcr versus t for A = 2 mm, where 
the approach velocity is varied by a factor of 17. In this 
case, the closed symbols represent v = 0.0001 A/ms, and 
the open symbols represent v = 0.0017 As/ms. The data 
for V = 0.0004 A/ms, shown previously, is within error 
of the V = 0.0001 A/ms data. 

Fig. 11 shows Cinit versus r for v = 0.0001 A/ms, v = 
0.0004 A/ms, and v = 0.0017 A/ms, where A = 2 mm. 
Three independent coalescence events are shown for each 
velocity. For v = 0.0001 A/ms and v = 0.0004 A/ms, 
the capacitance is very reproducible, and Cinit = 0.52 ± 
0.14 pF. However, at t; = 0.0017 A/ms, the capacitance 
is significantly larger and also less reproducible between 
different events. 

We do not see a significant change in the data when 
we increase the approach velocity up to t; = 0.002 A/ms 
for A = 1 mm and up to w = 0.0004 A/ms for A = 2 
mm. The behavior we see for drops A = 2 mm, v = 
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FIG. 10: Varying approach velocity, (a) Rcr versus r for 
A — 1 mm. The closed symbols show v — 0.0004 A/ms. The 
open symbols show v = 0.0020 A/ms. Data shown is the 
average of 12 individual coalescence events for each v. Three 
coalescence events occurred at each of four different /, from 
/ = 10 kHz to 10 MHz. (b) Rcr versus r for A = 2 mm. The 
closed symbols show v = 0.0001 A/ms. The open symbols 
show V = 0.0017 A/ms. 



0.0017 A/ms is consistent with increased deformation of 
the drops. A highly deformed drop would coalesce at 
the same c? as a less deformed drop, but tc would occur 
later, and the effective A could be expected to be larger. 
This leads to a lower i?cfl than a hemispherical drop for 
r > Tfiat, which we see for A = 2 mm and v — 0.0017 
A/ms when compared to v = 0.0001 A/ms. A large 
deformation would also increase Cinit significantly, which 
we also observe in the highest-velocity data for A = 2 
mm. The lack of reproducibility of Cinit for v = 0.0017 
A/ms may indicate that the radius of the flattening is 
not consistent between different events at large velocity. 



IV. CONCLUSIONS 

In conclusion, an electrical method has been used to 
study the coalescence of two salt water drops. This 
method allows us to observe an unexpected asymptotic 
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FIG. 11: Varying approach velocity, dnit versus r. Closed 
symbols show data from three separate coalescence events 
for V = 0.0004 j4/ms. Open symbols show data from three 
separate coalescence events for v = 0.0001 A/ms. The 
crosses show data from three separate coalescence events for 
V = 0.0017 A/ms. AH data shown was measured at / = 10 
MHz. 



regime which becomes visible at r < 10 ^s. Our data is 
consistent with the coalescence of two slightly flattened 
hemispherical drops. This is contrary to previous expec- 
tations, in which the drops were expected to maintain 
shapes described by quadratic minima. In addition, when 
A is varied by nearly a factor of three, we continue to ob- 
serve behavior which is consistent with the coalescence 
of two slightly flattened drops. Within error, we do not 
see a significant change in the data when we increase the 
approach velocity. 

A previous theoretical description [T3] has suggested 
that coalescence may occur as capillary waves cause re- 
peated connections of the gap between the two drops. 
Each connection would entrain a toroidal bubble of the 
outer fiuid. We see no evidence of this behavior in these 
experiments, which would appear as discrete jumps in 
Rcr during coalescence, as the neck width widens at each 
connection. It is possible that this behavior occurs on a 
timescale that is faster than the experiments described 
here are able to resolve. 

Previous experiments using high-speed imaging have 
been unable to resolve this early-time regime. Addition- 
ally, they found that electrical contact occurred 20 to 80 
/iS before the first motions of coalescence were observed 
visually [TSJ US]. If the drops are coalescing in a fiat region 
with d \QQ nm, it would be impossible to observe this 
stage of coalescence visually. It would only be possible to 
observe coalescence when d increased, entering into the 
hemispherical regime. We find that this occurs between 
1 and 10 /is after the initiation of coalescence, consistent 
with these observations. These previous experiments also 
postulated that coalescence occurred over a finite region 
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of radius ~ 100 /zm. In this case, Rcr would increase 
suddenly at the instant of coalescence, contrary to our 
observations. 

There are several possible reasons for the existence of 
a flattened region, and I suggest two here. One possibil- 
ity is that the flattening is an air effect. It has recently 
been shown that for drops splashing on dry surfaces, air 
plays a role in the dynamics of the impact [25' . Two 
drops approaching one another at finite velocity might 
trap a layer of air between them, which could have unex- 
pected consequences. A second possibility is the presence 
of surfactant. Although precautions were taken to avoid 
contamination [28], surfactants might still be present in 
small quantities. It has been seen that even a very small 
amount of surfactant can prevent a drop from coalescing 
with a flat fluid surface [25]. The repulsion due to this 
could explain the observed flattening. 

Understanding the flattening of two fluid drops as they 
approach each other could not only affect the many in- 
dustrial applications that rely on droplet coalescence, but 
also could illuminate other important physical questions. 
The origin of the thin film rupture which triggers coa- 
lescence is an active field of research j27j. The observable 
flattening of the drops could contribute to an understand- 
ing of this rupture. Also, as we begin to study the topo- 
logical changes that occur in microfluidics, behavior at 
the smallest scales and the earliest times is essential. 

The earliest stages of a topological transition in a fluid 
are when the analogy to a critical thermodynamic phase 
transition ought to be most accurate. As in drop break 
up, near the coalescence transition, the small-scale flows 
decouple from the large-scale flows. However, we observe 
that the bridge radius between the two drops scales differ- 
ently depending on the overall drop shape. The geometry 
of the system is crucial, a situation that does not have an 
analog in thermodynamic phase transitions. By studying 
these fluid shape transitions, we widen our understanding 
of the many unexpected ways in which nature produces 
these remarkable transformations. 

I am grateful to S. R. Nagel, L. N. Zou, X. Cheng, J. 
L. Wyman, N. Keim, E. Corwin, and J. Royer for helpful 
discussions and feedback. This research was supported 
by NSF MRSEC DMR-0213745 and NSF DMR-0652269. 



APPENDIX A: CHECKS ON METHOD 



Our experimental method, although similar to previ- 
ous methods, has many novel elements. We present here 
several checks that we have performed to validate our 
data. This includes isolating electrical effects as well as 
using our method to study a better-understood topolog- 
ical change, that of drop snap-off of water in air. 




FIG. 12: Varying electrical parameters in resistance measure- 
ments. Each set of data shown is the average of 12 individ- 
ual coalescence events, (a) Varying |Knai|. Rcr versus r is 
shown. Closed symbols show |Knai:| = 25 mV. Open symbols 
show |Vmaa:| = 50 mV. Crosses show |Vmaa;| = 250 mV. In- 
set shows Cinit versus r. Open symbols show |V;,iaa;| = 50 
mV. Crosses show IVmaa;! = 250 mV. (b) Varying DC com- 
ponent. Rcr versus r is shown. Closed symbols show no 
added DC. Open symbols show an added DC component of 
140 mV. Crosses show an added DC component of 35 mV. 
(c) Varying ionic concentration. Rcr/Pt versus r is shown. 
Closed symbols show a 26% by mass solution (saturation), 
o — 0.225 1/Q- cm. Open symbols show a 10% by mass solu- 
tion, a = 0.126 l/n- cm. 



10 



1. Varying electrical parameters 

a. Varying Voltage Magnitude 

We begin by verifying that our measurements do not 
change if the magnitude of the AC voltage across the cell 
is changed. For \V\ > 50 niV. wc cannot approximate 
the electrodes as a resistor in parallel with a capacitor, 
as the charge transfer reaction no longer gives us / ~ 
V. In addition, large \V\ could deform the shape of the 
drops by increasing the attraction or repulsion between 
the surfaces. 

During coalescence, as Zceii changes, the voltage across 
the cell, Vceii, also changes in magnitude. We can deter- 
mine a maximum | Vceii |, |Kriaa;|, by assuming that the 
full voltage drop supplied by the function generator oc- 
curs across the experimental cell. We vary IKnaKl from 
12.5 mV to 500 mV, and find no significant difference in 
our data, as shown in Fig. 12(a) for three sample |Knax|- 
The noise becomes comparable to our signal at 12.5 mV, 
and all data used in the analysis is taken at |V^ax| = 50 
mV. 

Additionally, as shown in the inset to Fig. 12(a), we 

examine the effect of varying iVmaa;! from 50 mV to 500 
mV on the measurement of Cinu using measurement fre- 
quencies from / = 6 MHz to 15 MHz. No difference is 
seen outside of error. 



6. Varying DC component of input signal 

In addition to the AC signal applied across the exper- 
imental cell, a small amount of DC signal is observed. 
Any DC signal applied across the experimental cell be- 
fore coalescence will polarize the cell, effectively charging 
it up like a capacitor. At the instant of coalescence, this 
capacitor discharges, and a DC spike is observed. This 
spike has a typical maximum size ^ 1 — 10/iV, which is 
less than 1% of the typical output signal. We average 
the signal and remove this DC contribution before anal- 
ysis. We check the validity of this by explicitly adding 
up to 140 mV of DC to our signal. This did not alter our 
measurement within error, as can be seen in Fig. 12(b). 



c. Varying ionic concentration 

All data presented in the main body of the paper was 
taken with NaCl in water at saturation, or 26% by mass. 
We took the same set of measurements with solution of 
NaCl 15% by mass and 10% by mass. The conductivity 
a changes by a factor of two between the NaCl solution 
at saturation and that at 10% by mass. All resistances 
measured should be inversely proportional to a, and so 
to compare the different solutions we looked at Rcr<7- 
RcrP' versus t is plotted in Fig. 12(c); and does not 
change, within error, as the concentration is varied. 



With changing concentration, the fluid parameters also 
change slightly. For a solution of 10 % NaCl by mass, p = 
1.0707 g/cvc?, V = 1.115 cSt, and 7 = 76.05 dyne/cm. 
Calculating the predicted prefactors for the model, we 
find that the differences due to the fluid parameters are 
minimal, and we should not see effects from them outside 
of our experimental error. 



2. Drop snap-off 

Finally, we used our method to study drop break-up. 
We compare the output of our method against previous 
work, as well as our own calculations, and find that they 
are consistent. 

It has been previously observed that during drop snap- 
off, the neck between the two drops forms a self-similar 
cone with an angle of approximately 18°. We numeri- 
cally calculated the resistance of a truncated cone of fixed 
larger radius as a function of the smaller radius, r„ec*;, 
and found Rcr = 1. 18 /ar neck where a is in units oil/il- 
cm. 

We used a Phantom V.7 fast digital camera at 144,000 
frames per second to measure r„ecfe as a function oi t — 
to = T, where to is the instant of snap off. As shown in 
Fig. 13(a), the best fit to our data with r < 1 ms yields 
rneck = (1.2±0.6)r-°-'^^=^°-°<^, which is in agreement with 
previous measurements of drop break-up of water in air. 

As shown in Fig. 13(b), Rcr versus t for a snap-off 
event yields a best fit of Rcr = (5.7 ± 0.9)r-°-^^=t°-°i. 
This is shown as the dashed line in the figure. Combin- 
ing the electrostatic calculation with the measured de- 
pendence of Vneck on T. wc predict that Rcr — (5.24 ± 
2.63)t-° *5^±" "^', within error of the data. 

This indicates that the scaling law for break up of wa- 
ter in air persists to timescales of 10 ns. Using the rela- 
tion above, we calculate rneck from Rcr- rneck versus t 
is shown in Fig. 13(c). The solid line is the best fit to 
the data, r„ecfe = (0.75 ± 0.01)t-°-6^±° °^ 

APPENDIX B: CAPACITANCE 
MEASUREMENTS 

Before the bridge is formed, we need to separate the 
capacitance of the coalescing region, Ccr from the capac- 
itance of the total arrangement of conductors. As stated 
in the experimental section, in order to achieve this, I 
measure the capacitance of the cell before any drops are 
formed on the nozzle tips, Cceih and subtract this from 
the Ccr measured in the last /xs before coalescence oc- 
curs. 

I justify this approximation by simulating the axially 
symmetric part of the cell, including the nozzles and the 
tubes filled with NaCl solution. I find that for a fixed 
separation of the two drops (^ = 1 mm) of d = 1.5 mm, 
Ccr — 0.45 pF from the simulation. Simulating the cell 
with no drops formed gives Cceii = 0.44 pF. 
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For a 1.5 mm separation between the drop tips, the 
measurements show Cqr = 0.66 ± 0.02 pF, and Cceii = 
0.62 ± 0.02 pF. When there is no fluid in the cell, we 
measure a stray capacitance of 0.2 pF, which accounts for 
the difi^erence between the measurement and simulation 
if it can be considered to be in parallel with Ccr- 



The capacitance of the cell outside the "coalescing 
region" is in parallel with Cinit, thus, I estimate that 
Cinit = Ccr — Cceii- The simulation when d = 1.5 mm 
yields Ccr — CceH = 0.01 pF, which is consistent with 
our measurement of Ccr — Cceii = 0.04 ± 0.03 pF. 
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FIG. 13: Drop snap off. (a) Tnecfc versus r. r„ecfc was mea- 
sured using high-speed imaging at 144,000 frames per second. 
Tiie solid line shows r = 1.2t'"'^"' (b) Rcr^ versus r. The 
dashed line shows Rcr = 5.7t~'^'''^ (c) r„eck versus t. rneck 
calculated from electrical measurement of Rcr- The solid line 
shows rneck = 0.75t-°-^'' . 



In order to find Cinit just before coalescence, I measure 
Cceii at the nozzle separation used for the coalescence 
measurements, and subtract this from Ccr- For A = 1 
mm and a nozzle separation of 2 mm, Cceii = 0.89 pF. 
The simulation predicts Cceii = 0.64 pF in this case. This 
measurement is repeated for A = 2 mm and d = 4 mm, 
where Cceii = 0.92 pF, and at ^ = 0.75 mm and d = 1.5 
mm, where Cceii = 0.81 pF. 
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